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Abstract 

We discuss a few possible strategies for measuring the polarization of the 
Ab baryons produced in e"'"e~ -annihilation at the Z resonance through their 
inclusive semileptonic decays. After reviewing the existing methods, an ex- 
tension is proposed, based on the ratio of the averages of the squared electron 
and neutrino energy, including both perturbative and nonperturbative cor- 
rections. This variable minimizes the statistical error on the Ab polarization, 
while keeping the systematic theoretical errors at the level of 1-2 %. A number 
of other polarization-sensitive variables are also discussed, such as averages of 
ratios of the electron and neutrino energy and the distribution in the difference 
of the electron and neutrino rapidities. 
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I. INTRODUCTION 



It is a well-known fact that the bottom quarks produced in e^e~-annihilation at the Z 
peak are predicted to be very strongly polarized. This can be understood by noting that at 
this energy the production process is practically dominated by Z-exchange, whose coupling 
to quarks is proportional to T3 — sin^ Oy/Q. In the limiting case when 9\y = 0, the right- 
handed quarks (for which T3 = 0) decouple completely from the neutral current and thus 
the produced quarks are purely left-handed. Increasing the value of the weak angle 6w to 
its physical value somewhat diminishes the effect (more for the u-type than for the d-type 
quarks), such that for sin^ 6'vi/=0.23, a polarization degree P = —0.936 is expected for the 
d-type quarks. 

The b-quark mass effects are suppressed by a factor of m^/m| and are very small, less 
than one tenth of a percent. The radiative corrections to this prediction have been computed 
to one-loop order and have been found to decrease P slightly, by about 2%. Also, the 
polarization depends only weakly on the kinematical details of the production process. 

Eventually the b-quark is observed through its hadronization products. The b quark 
hadronization into mesons (even excited ones) deletes any memory of the original quark 
polarization . On the other hand, if the quark ends up as a Ab baryon (in which the light 
constituents combine to spin 0), the latter is expected to retain a large part of the initial 
polarization of the quark [^. 

A measurement of the polarization of the produced Ab's could therefore help to test 
directly this prediction of the Standard Model and obtain information about the details of 
the hadronization process. A first measurement is already available 0, which gave an in- 
triguing ly small value P = -0.23lg:|. This measurement made use of the method proposed 
in [0, based on the ratio of the average electron and neutrino energies y = {Eg) / {E^) in the 
laboratory frame. This is just one of the many approaches which have been given which are 
sensitive to Ab polarization effects, based on its semileptonic decays to charmed hadrons. 
Thus, in the average lepton energy in the laboratory frame has been suggested as a 
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measure of the Ab polarization. Unfortunately, this quantity is rather sensitive to the exact 
values of the quark masses, which are known only imprecisely. 

In contrast, the variable y = (Ee) / {E^) proposed in ||^ presents the advantages that the 
poorly known b-quark mass drops out and the dependence on the ratio m^/m^ is compen- 
sated to a large extent. Also, the dependence on the details of the fragmentation process 
(the fragmentation function) disappears almost completely, a problem which had been cir- 
cumvented in HJ^ by a combined use of low-energy experimental data. Recently, yet a 
different method has been put forward in [0. 

The purpose of the present paper is twofold: a) to give a detailed analysis of the theo- 
retical errors which affect an extraction from experiment of the Ab polarization and b) to 
discuss a number of other variables which can be used to measure the polarization. We find 
that the theoretical uncertainty is dominated by the one in the ratio of the quark masses 
p = m1/m^ and by the uncertainty in the strong coupling as- In section 2 we discuss a 
generalization of the variable proposed in |]^, the ratio of the moments of the electron and 
neutrino spectrum. We argue that the ratio of the second moments minimizes the statistical 
errors in the polarization, while being also insensitive to theoretical uncertainties. In sec- 
tion 3 two other polarization-dependent variables are proposed, the averages of ratios of the 
electron and neutrino energies and the decay distribution in the difference of the electron 
and neutrino rapidities. 

II. MOMENTS OF THE LEPTON SPECTRA 

The polarization of a Ab baryon affects directly the spectra of the leptons (e and v) 
produced in its semileptonic decays. The qualitative nature of the change can be obtained 
by examining the angular distribution of the leptons in the rest frame of the decay: the 
electrons tend to be emitted antiparallel to the spin of the Ab, whereas the z/'s are emitted 
preferentially parallel to the spin. In the laboratory system this translates into a larger 
number of electrons being emitted coUinear with the jet containing the decaying b quark 
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(assuming negative polarization) and hence more energetic as compared with the unpolarized 
case (and the opposite for neutrinos). Consequently, the electron spectrum in the laboratory 
frame will be hardened in the presence of negative Ab polarization, contrary to the energy 
spectrum of the neutrinos, which is softened. 

It is more convenient to study the moments of the lepton spectra, rather than the spectra 
themselves. The theory of the inclusive semileptonic decays of heavy hadrons predicts that 
the lepton spectra are given, in a first approximation, by the free-quark decay results. The 
corrections to these results are of order l/m^ [ITTHTsU and can be parametrized in terms of a 



few matrix elements. We write the general form of the electron spectrum in the rest frame 
of the decay as 

dr \ v P 

dxedcos^e = ^!(^^l(F^(^^^P) + P^-'(^M^^^P)) ' (2-1) 

where = E2/2mi, is the reduced electron energy, 9^ is the angle between the electron 
momentum and the Ab momentum, p = m'^/m'^, and P is the polarization. In writing this 
we have assumed that the Ab is longitudinally polarized, which is true to a high degree of 
accuracy^. Thus P is antiparallel to the Ab momentum. 
To leading order in 1/mb one has 

Feoixe, p) = ""^^^"^"^'^' [(1 - Xe)(3 - 2Xe) + p(3 - x^)] (2.2) 

Xe) 

MXe, P) = ^^^J— ^^|^[(1 - X,){1 - 2X,) - p(l + Xe)] . (2.3) 

The l/m^ nonperturbative corrections to these relations have been computed in |T2| , p!3| and 
the a^-order correction in [0,0] . First we will neglect all corrections and use the free-quark 



decay results. The energy of the electron in the laboratory frame is given by a boost 

e, = ^{e; + PpI) (2.4) 



^The transverse polarization, averaged over the orientation angle with respect to the e~^e beam, 
is predicted to be less than 1% in the Standard Model |14,15|. 
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where center-of-mass quantities are starred. With the help of this relation the moments of 
the laboratory frame electron energy can be reduced to the moments of the two functions 

r-l-p 



1"P 



(2.5) 
(2.6) 



n 




K) 

(7") 
\ / / 


1 


7 - P 19 - 7P 
ID ^ ID 





2 


8(4 -P) 16(11 -5P) 
45 ^ 45 




3 


2(3 -P) 2(25 -13P) 
7 P 7 


4(5 + 3P),, 
35 


4 


8(5 - 2P) 64(7 - 4P) 
35 ^ 35 


«(2 + 2P)(^_3) 
35 ^ 


5 


4(11 - 5P) 20(31 - 19P) 
27 27 


252 ^ 


6 


128(2 -P) ^256(17 -IIP) 
105 ^ 105 


105 ^ 



Table 1. Average values of the first six moments of the electron (second column) 
and neutrino energy (third column) in the laboratory frame, to first order in p. 



To first order in p they are given by 
„ 6 + n 



J p. 



(3 + n)(4 + n) 
2 + n 



3p 



+ 3p 



3 + n 
2 + n 



(3 + n)(4 + n) ' 3 + n 
5 



+ 0{p' 



(2.7) 
(2.8) 



In Table 1 the moments of the reduced lepton energy in the laboratory frame are listed 
as functions of P, to first order in p. {(3) has been taken equal to 1 in computing these 
expressions, since the b-quarks produced at the Z pole are practically ultrarelativistic. The 
qualitative behaviour of the P-dependence is just what one expects from the general con- 
siderations mentioned at the beginning of this section: the moments of the Xg-distribution 
grow as P approaches the completely polarized value —1. The relative increase with respect 
to the unpolarized case grows with n, due to the fact that the higher moments increasingly 
probe the region of high values of Xe, where the spectrum is hardened in the polarized case. 

Unfortunately, the higher moments display also a larger sensitivity to p, which is the 
most important source of uncertainty in these calculations. Previous works used a rather 
broad range of values 0.06-0.12 @J^J3I f*-*^ "which obtains when the individual quark masses 
are varied independently within the limits = 4.8 ±0.3 GeV and = 1.35 ±0.15 GeV. It 
is worth noting that a simple use of the heavy mass expansion in the HQET can be used to 
narrow down this interval. This can be done by writing simultaneous heavy mass expansions 
for the masses of the corresponding heavy mesons 
1 «2 

-{mB + 3mB*) = mb + A + i^ + 0{l/ml) (2.9) 
4 2mb 

1 

-{ruD + SniD*) =m, + A + ^ + 0{l/ml) . (2.10) 
The binding A and the heavy quark kinetic energy p^ are still imprecisely known parameters. 



For example, in the bound A > 500 MeV was obtained and in [|T8| correlated bounds 
on A and /i^ are given, as well as the lower bound A > 410 MeV. Lattice calculations give 



much smaller results pO[ for A. We will use as a conservative number A = 450 ± 50 MeV. 



A QCD sum rule calculation gave the value p^ = (0.60 ± 0.10) GeV^, in agreement 



with model-independent bounds derived in |23,24|. We will use this central value but double 



the error bars. On the left-hand side of ( p.9|J2.10| ) we use average masses over the respective 



isospin doublets, which gives liniD + 3mD*) = 1973 MeV and i(mB ± 3mB*) = 5279 
MeV. Neglecting the terms of order 1/m^ and higher in ( p^.9[ - f^.l0] ) one obtains (mb, rric) = 
(4.745 GeV, 1.355 GeV) ^(4.871 GeV, 1.547 GeV) and correspondingly p = 0.0815-^0.1009. 



These numbers are similar to those quoted in i.e. mb = 4.80 ± 0.03 GeV and rric = 
1.35 ± 0.05 GeV (see also the references cited in [p^). 

Returning to the analysis of the p-dependence of the moments of the electron energy 
spectrum, we present in Table 2 the values of these moments corresponding to the fully 
polarized Ab case (P = —1) and in the unpolarized case (P = 0). In computing these 
numbers the dependence on p has been fully taken into account, not only to first order as 
in Table 1. 



n 


(7") 


{<) 


P = -1 P = 


P = -1 P = 


1 


0.648;°:°- 0.583;°:°!° 


n q4fi-0-005 n c;iQ-0.007 
U.O^U_^0.005 '-'■'-'-'-'^+0.007 


2 


5q2-°-°2^ 4qfi-o oi6 


200"°-°°^ 400-° °^^ 


3 


632-°-°^3 f^ns-°-°2'^ 


Q -I4Q-0.006 r, 070-0.015 
U.i^J_l_0,oo6 '-'■'-"'-'+0.016 


4 


0.741;°:°^° o.566;°:°i 


n 1 on-0.007 n oqi -0.021 
U.±OU_^0.007 '-'■'-''^J^ +0.022 


5 


0.924;°:°^^ 0.684;°:°^? 


Q . 27-0-008 Q 444-0.029 
'-'•-'-'^'+0.009 U-^^^+O.OSl 


6 


-1 9r,7-0.115 r, 070-0-078 
l.ZU(_,_o.i29 U.O(U_^Q 087 


M 1 qq-O.OlO M c-qq-0.041 
U.lOO_|_0.011 "-'•^'-''-'+0.045 



Table 2. Sensitivity to P versus sensitivity to p in the electron and neutrino 
spectrum. The value p = 0.091 ± 0.010 has been used (see the explanation in 
the text). 

The result is that, for the first moments, the error due to p is only marginally smaller 
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than the sensitivity to P With growing n both of them increase such that for larger 

n the uncertainty induced by p is offset by the increase in the sensitivity to P and the 
efficiency of the variable becomes better. However, the additional uncertainty connected 
with the fragmentation-dependent factor (7"), the need of considering high moments in 
which the statistical errors become appreciable and the comparatively small efficiency make 
the use of the electron spectrum little attractive as a polarization analyzer. 

Similar considerations can be made for the case of the inclusive neutrino energy spectrum. 
This is described by a formula analogous to with the corresponding functions F^, and 
Ji/ reducing to the well-known free-quark decay expressions in the infinite mass limit 

Fu{x,,p) = ■Mx,,p) = ^^K^ - P - ^'^)\ (2.11) 

1 - Xi, 

Their moments can be expanded to first order in p with the result 

Fun = J.n = , - ^ + 0{p') . (2.12) 

(?2 + 3)(n + 4) n + 3 

The results for the first few moments of the neutrino reduced energy in the laboratory frame 
are shown in the third column of Table 1. Their relative sensitivity to P and p are given 
in the last two columns of Table 2. One can see that the main conclusions drawn in the 
electron spectrum case can be extended to the neutrino case. 

More satisfactory variables are obtained by taking the ratios of the moments of the 
electron and neutrino spectrum. Numerically, they are much more stable against changes 
in p because the electron and neutrino moments go in the same direction under a variation 
in p. At the same time, the sensitivity to P is strongly enhanced (see Table 3). Most 
importantly, the boost-dependent factors (7") cancel out in the ratios i/n, which are therefore 
independent of details of the fragmentation process [0. 

The sensitivity to P increases rapidly with n, from a relative change of 60% in yi as P 
goes from to 1 to over 450% for y^. In [0] this relative change has been used as a measure 
for the sensitivity of a given variable to the polarization. It is clear, however, that the effect 
of the statistical errors on each particular variable has to be taken into account as well. 

8 



We propose instead as a more realistic measure of the relative efficiency of two different 
variables, the ratio of the errors on P induced by each of them with identical statistical 
data. 



n 


- {xn) 


SPi 


P = -l P = 


P= -1 P = 


1 


1 Q7Q-0.009 -1 -19/1-0.003 
I.O/O^QQIO -'-■-'-'^^+0.003 


1.000 1.000 


2 


9 QC^7-0-026 1 900-0.007 
z.y<J/ +0.028 -'-■^'^^+0.007 


0.81472 1.05572 


3 


A OQQ-0.051 1 o/i 7-0.010 


0.89373 1.35273 


4 


rr «C4-0.085 447-0.014 
O.UO^_,_0.o90 -'-•'±'±'+0.015 


1.03774 1.75674 


5 


7 904-0127 P-49-O.OI8 
'■^^^+0.135 ±•'-'^'^+0.019 


1.21975 2.25375 


6 


Q n'^'^-O l''"'' 1 fiQ9-0.022 


1.43176 2.84276 



Table 3. Ratios of the moments of the electron and neutrino energy spectra. 

Explicitly, assuming that the variable y„ can be measured with a statistical error Si/n, 
this allows for a determination of P with an error 



dp^ 



(2.13) 



The error of y„ is, in turn, given by 

ny, 



/p2{n-l)v ,p2(n-l)v 
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(2.14) 



where is the number of events and 5Ee,u are the statistical errors with which the electron 
and neutrino energies are measured. Assuming that the systematic errors dominate over the 
statistical ones, we will take 6Ee ~ 6E^. Using the data in Table 2 for the average energies 
one obtains the numbers shown in the last two columns of Table 3 for the ratio of the errors 
in P induced by yn and yi. The boost-dependent factors 7„ are given by 

^" = ^^\^ (2.15) 

and are of the order of the unity. For example, by using a Peterson ansatz for the fragmen- 
tation function of the form f{z) = with e = 0.016 the following average values 
are obtained: (7) = 6.244, 72 = 0.981, 73 = 0.989 and 74 = 1.009. 

One can see that, as far as the statistical errors are concerned, y2 allows an extraction 
of P with an error at least 20% smaller than yi (note that 72 < 1 independent of the 
fragmentation model adopted) in the fully polarized case (P = — 1). In the vicinity of P = 
the two variables are equally sensitive to statistical errors. 

Incorporating the nonperturbative corrections amounts to replacing the free-quark decay 



moments (Ora according to Pi ^ Pi, Ji ^ (1 + eb)Ji, P2 ^ (1 + iK^)F2, J; 



2 



(1 + eb + |P'b)^2 ||T3[- The nonperturbative matrix elements appearing in these relations 
have the form P'b = fil/2m1 and eb = fil/ml. We will use the value ^"j. = — suggested 



by a model-independent bound |2^ and fi'^ = 0.6 ± 0.2 GeV^ as already explained above. 

The analytic expressions for the radiative corrections can be found in 0] . We have com- 
puted the corresponding corrections to the moments of the electron and neutrino spectrum 
numerically. The value of the strong coupling has been varied between 0.1 and 0.3. The 
results, including the errors associated with the variation of the parameters involved, are 
presented in Table 4. 
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variable 


central 


p 

r 




as 


Total 




1.131 


+0.003 
-0.003 


+0.000 
-0.000 


+0.004 
-0.003 


+0.007 
-0.007 




1.901 


+0.011 
-0.010 


+0.004 
-0.004 


+0.011 
-0.010 


+0.027 
-0.023 






+0.005 
-0.005 


+0.003 
-0.003 


+0.004 
-0.003 


+0.013 
-0.011 


Vo( P = 0) 


1.253 


+0.007 
-0.007 


+0.000 
-0.000 


+0.007 
-0.008 


+0.016 
-0.014 


VoiP = -1) 


3.039 


+0.039 
—0.029 


+0.011 
-0.011 


+0.031 
-0.028 


+0.076 
-0.065 




2.425 


+0.011 
-0.010 


+0.009 
-0.009 


+0.009 
-0.008 


+0.030 
-0.026 



Table 4. The contributions of the different sources of theoretical errors to yi 
and I/2- Ry„{P) = yn{.P) / VniS^) are the ratios ?/„ normalized at 1 at P = 0. 

The systematic theoretical errors can be further reduced by considering the ratios 
Pyn{P) = yniP)/yniP = 0) p. Since these ratios differ from ?/„ only by a constant, the 
arguments following after Eq.(15) concerning the statistical errors apply to them without 
changes. Assuming zero statistical errors, the total theoretical uncertainty of -Ry^ is reflected 
into an error of 2.3% in P at P = —1, compared with about 4.3% when considering yi at the 
same point. In a similar way, Ry^ induces an uncertainty in P of 2.0% at P = —1, compared 
to an error decreasing from 3.9% at the same point to 3.3% at P = for ?/2- In the Figure |l|, 
the Ry^ and Ry^ variables are plotted as a function of the Ab polarization for p = 0.09. 
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FIGURES 




0. -0.25 -0.50 -0.75 -1. 

PAb 



FIG. 1. The Ry-^ and Ry^ variables versus the At, polarization for p = 0.09. 

In conclusion, y2 and Ry^ could prove to be useful as variables for measuring the Ab 
polarization in addition to yi (respectively Ry^), due to their reduced sensitivity to statistical 
errors. 

III. OTHER POLARIZATION-DEPENDENT VARIABLES 

In this section we discuss a few other variables which can be used to measure the polar- 
ization of the Ab baryons. For example, one could consider the average ratio of the electron 
and neutrino energies. In terms of rest-frame quantities this can be written as 

We will keep (3 fixed for the moment and only average over it at the very end. Let us 
compute the average ( |3.1| ) in the free-quark decay approximation. The differential decay 
rate dF is given by 

12 



1 2 



XgXj^(l + P COS 6i,)dxyd cos O^d cos Oeydcj) . 



(3.2) 



where 



I- p-x^ 



1 + f (cOS^e 



(3.3) 



and /(p) = 1 — 8p + 8p'^ — P^ — 12p^ log p. is the angle between the electron and neutrino 
moments in the rest frame of the decay and denotes the angle between the decay plane 
and the (pj,, p^) plane in the same reference frame. After integrating over all variables in 
( p.l|) the following result is obtained 



Vd 



6 



/(P) 
1 



{-^(1 - P)(l - lip - 47p2 - 3p') + p'(3 + 2p) logpj 



(3.4) 



_ p (3-P l + (3' 



;i - p)(l - 8p - 17p2) - 6p2(3 + p) logp 




0. -0.25 -0.50 -0.75 -1. 



Ab 



FIG. 2. Ry^ and Ry^ as functions of the Ab polarization. 
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Since (3 is very close to 1, it can be replaced everywhere with this value, except in the 
argument of the logarithm, which can be written as 

logi±| = 21og27 + 0(/5-l). (3.5) 

With this approximation, (|3.4|) becomes 

VD = {-^(1 - P)(l - lip - 47p2 - 3p3) + p\3 + 2p) logpj (3.6) 



+ T7T + (1 - ^)(log27)) Rl - p)(l - 8p - 17p2) - Qp\3 + p) logp 



2/(p) I 12' 
2 

7(p) 

The fragment at ion- dependent quantity (log 27) can be eliminated by considering also 
the average of the inverse ratio 



This can be computed in a similar way to yn with the somewhat lengthy result 
I// = {-^(1 - P)(3 - 29p - 65p2 + 31p3) + p2(6 - p^) logp} 
+ ^^jf^ [1(1 - P)(2 - 13p - 4p^ + 3p3) - 2p^(3 - p) logp} 
+ {-^(1 - P)'(l + lOp + p') - p(l - p') logp} (3.8) 

+ ^7^^ {^(1 - - 8P - 17p^) - 2p^(3 + p) logp} 
- {^(1 - P)(l - 5p + 13p^ + 3p3) + 2p3 logp} . 

To a good approximation, the last term in (|3.8| ) can be neglected because of the addi- 
tional suppression by a factor of 1/7^. An estimate performed with the help of the Pe- 
terson fragmentation function used in the preceding section gives (log 27) ~ 2.497 and 
((log 27)772) ~ 0.064. 

The variable yo should be used carefully; due to detector smearing the neutrino energy 
can be close to zero and create overflows. The solution to this problem is to do a cut 
on neutrino energy around zero. This cut together with other analysis cuts {e.g. lepton 
identification cuts on the lepton momentum at 3GeV and on its transverse momentum at 
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IGeV/c) will bias the values of yo and yi. The normalizations Ryj^ = yo^P) /yoiS^) and 
Ryi = yi{P)/yi{^) can be used to eliminate the acceptance and reconstruction effects (see 
Figure 0) . A Monte Carlo study shows that the relative errors obtained in a sample of N 
events are the following 



Pa. = 0. Pa 



b 



cTR,^ 2.4 2.0 



Ryo VN VN 
any, 2.4 2.8 



Ryi VN VN 



A simultaneous fit to ( p.6|) and (|3.8| ) will give the polarization P. The fragmentation 
variable (log 27) can be eliminated between Ry^ and Ry^ and the polarization value extracted 
as a function of both neglecting the term proportional to (^^^1^). 



The fragmentation-dependence of the variables yD and yj (pX,3.8) can be traced back 
to the fact that the energy ratios in and respectively ( |3.7| ) are themselves not Lorentz- 
invariant under boosts along the direction of motion of the Ab. One can construct boost- 
invariant distributions and averages (hence insensitive to the details of the fragmenta- 
tion process) by considering quantities which are explicitly invariant under such Lorentz- 
transformations. 

One such quantity, which turns out to be also sensitive to Ab polarization, is the difference 
of the rapidities of the electron and of the neutrino along the boost direction 

A = Ce-C.- (3.9) 

In general the rapidity of a particle is defined as C = | log with E its energy and q 
the longitudinal momentum. For a massless particle the rapidity can be also expressed 
through the angle 6 between its direction of motion and the boost axis, tanh(^ = cos^. We 
recall that in the polarized case (P = —1), the electrons tend to be emitted preferentially 
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along the boost direction and the neutrinos in the opposite direction. This means that the 
distribution in A should be shifted towards positive values for a polarized Ab compared with 
the unpolarized case (see Figure 




FIG. 3. The difference between cliarged lepton and neutrino rapidities A for two values of the 
Ab polarization P\^^ = — 1. and P\^^ = 0. 

The simplest way to obtain the decay distribution differential in A is to parametrize the 
final state of the system from the very beginning in terms of the set of coordinates {(, r, 0) 
for each particle: the rapidity (, the transverse momentum to the boost axis r and the 
azimuthal angle 0. This gives 

dF 



2{2n)^ ^e("^b - - 2mbre coshCe)^ 



X 



ml cosh(Ce + C^) + ml cosh(Ce - C^) + [ml - ml) cos( 



(3.10) 



X (cosh Cu + P sinh (^^ 



{mbCOshCi. + re[cos(0e - 0^) - COsh(Ce - C;^)]}^ 

The electron transverse momentum takes values from to mb(l — p)/(2coshCe)- It is 
therefore helpful to introduce a reduced transverse electron momentum r defined as r = 
2reCoshCe/'^b which will take values between and 1 — p. 
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Theoretical prediction 
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FIG. 4. (A) variable as a function of the Ab polarization simulated for five different values 
= 0.0, —0.25, —0.5, —0.75 and — 1.0 (points) and the theoretical prediction (continous line). 
The Monte Carlo samples contain 3000 events and the quark mass ratio has been taken p = 0.09. 



After integrating over 0e,i^ and transforming from (Ce, (^u) to C+ = Ce + Cu ^ i W^ 
the following result is obtained 



dr 



8(277)3 



r(l — p — r)^ < (1 — p) 



+ 2(S - r cosh A)^ 



1 - p - r)S 



[S - r cosh A] [Sr^ + 2(S - r cosh A)^ 



X [cosh (^_|_ + cosh A + P (sinh (^_|_ — sinh A)] drdAd^ 



+ ) 



(3.11) 



where S = cosh + cosh A and R = sinh^ A — 2rS cosh A + S^. 



The exact result obtained after integrating over r and C+ is very complicated, so we only 
give in analytic form the simpler expressions corresponding to the case of a massless final 
quark p = 0. In this limit the decay distribution has the form 



"^'-m^(/(A)-PsinhA^(A)) 



dA 4(27r)= 



(3.12) 



17 



where both /(A) and g{A.) are even functions of their argument. For positive values of A 
they are given by 

36s^ + 180s2 + 175 , A ^, A\ As^ + 218s^ + 715s^ + 525 



/(A) = -7(A)c 



+ (log(l + e^) - I) 



34s6 + (827 - 132c)s^ + 10(232 - 127c)s2 + 1575(1 - c) 



(3.13) 



and 



g{A) = -/(A) 
+ 



A\ 45^ + ISOs^ + 315 
+ log 1 + - - c 



32s^ + 130s^ + 105 

2^8 ' V'^^v- ' - / 3s8 

4(11 - 3c)s^ + 5(64 - 43c)s2 + 315(1 - c) 
6^8 ■ 



(3.14) 



We have denoted in these expressions 

1 



1(A) = — — f-2Li2(-e^) --A^- — 
^ ' sinhA I ^ ^ 2 6 



(3.15) 



and c = cosh A, s = sinh A. 



p 


c. 




0.000 


0.667 


4.089 


0.081 


0.633 


3.695 


0.091 


0.630 


3.662 


0.101 


0.627 


3.631 



Table 5. Average values of the difference of the electron and neutrino rapidities 
A = — Ci/, and of its cube A^, at P = — 1, for a few values of p. 

The singularity of /(A) and g{A) at A = is only apparent; at this point these functions 
have power series expansions of the form 

62 8 , .\ . ^ / 91 32 



/(A) = - Aiog2) - A^ {-^ + ^log2) 
' V945 315 ^ y V 2970 495 ^ ) 



-A' 



( 285209 17632 



V3243240 135135 



log 2 +C(A^ 



/AN A 73 32 , ^\ ^ 2 / 655 304 , ^ 
g{A)= (_-_log2) -A^l^-— log2 



V315 105 



2079 693 



V 117936 2457 ^ ^ > 



(3.16) 



(3.17) 
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As a measure of the polarization one can take the average value of an odd power of A. 
These are proportional to P and have the form 

(A") = -C„P (3.18) 

with Cn positive numbers. The first two coefficients C„ {n = 1,3), obtained by numerically 
integrating ( p.ll| ), are shown in Table 5 for a few values of p. In Figure ^ a Monte Carlo 



simulation was used to estimate (A) as a function of polarization. 

With the given error on p, (A) allows an extraction of the polarization with a theoretical 
uncertainty of the order of 2.1 % in the fully polarized case (and proportional to P for 
\P\ < 1), whereas (A^) gives a slightly larger error of 2.3 %. The statistical errors on the 
Ab polarization value are estimated in a Monte Carlo simulation to be 0.9/V^ with the 
number of events. While less precise than the variables discussed in section 2, nevertheless 
we believe that the methods presented in this section might prove to be useful as well, even 
if for cross-checking purposes only. 
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APPENDIX: 



We present here, for completeness, a few details concerning the radiative corrections to 
the first moments of the electron and neutrino spectrum. The corrections to the moments 
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of the functions Fe,u{,x) and Je,u{x) (describing the rest-frame decay distributions) have the 
form 

Fen -^en^ ^ / -^e fl ("^e) ; Jen ^en^ ^ / d^e-^eJl (-^e) 5 

TT Jo vr JO 

Fun J^i/n^ ^ / dXjy -^ufx (j'l/) ) Jun ^un^ ^ / X^j^ (•^i') • 

TT Jo TT Jo 

The polarization-dependent variables y and ?/2 are written in terms of these moments as 

3Fei -|- Jel-P 2Fe2 "t" Je2-P 



3-Fi/i -|- JuiP 



1/2 



2-F,y2 + Ju2P 



(Al) 



p 


^eO 


fo 


0.0815 


0.2760 


0.1173 


0.0912 


0.2577 


0.1074 


0.1009 


0.2406 


0.0984 



Table Al. The zeroth moment of the function Fe{xe) appearing in the unpo- 
larized 6-quark decay rate and its associated radiative correction. 



The functions fi (x) and ji (x) are defined in Eq.(16) and thereafter of Ref. [jT7 



p 


^el 


/r 


j(FQD) 
^el 


jT 




ft 


j{FQD) 


i+ 
Ji 


0.0815 


0.1638 


0.0701 


-0.0556 


-0.0152 


0.1453 


0.0651 


0.1453 


0.0659 


0.0912 


0.1503 


0.0630 


-0.0498 


-0.0137 


0.1337 


0.0588 


0.1337 


0.0595 


0.1009 


0.1380 


0.0566 


-0.0447 


-0.0124 


0.1231 


0.0531 


0.1231 


0.0537 



Table A2. The n = 1 moments and their radiative corrections. 

In Table Al we show the n = moment of Fe{xe) and its radiative correction, which are 
needed for the normalization of the absolute values of the higher moments. However, they 
cancel when considering the ratios of moments ?/„. 

The values of the free-quark decay moments F^^^^^ and J!^^^^^ and the moments of 
fi{x), if{x) (denoted as and j^) for n = 1 and n = 2 are listed in the tables A2 and 
A3. 
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p 




f2 


jiFQD) 
■Je2 


J2 


UFQD) 
^p2 


ft 


j{FQD) 
■Jv2 


Jt 


0.0815 


0.1062 


0.0461 


-0.0420 


-0.0121 


0.0852 


0.0403 


0.0852 


0.0407 


0.0912 


0.0960 


0.0407 


-0.0370 


-0.0106 


0.0774 


0.0359 


0.0774 


0.0362 


0.1009 


0.0866 


0.0359 


-0.0327 


-0.0094 


0.0702 


0.0320 


0.0702 


0.0323 



Table A3. The n — 2 moments and their radiative corrections. 

Prom these tables one can see that the radiative corrections change the n — 1 moments 
by 10-17% (for ag = 0.3). The corresponding effect in the n = 2 moments is somewhat 
larger, of 11-18%. 



21 



REFERENCES 

[1] J.G.Korner, A.Pilaftsis and M.Tung, Z.Phys.CGS 575 (1994). 
[2] M.Tung, Phys.Rev.'D52 1353 (1995). 

[3] S.Groote, J.G.Kdrner and M.Tung, MZ-TH-95-09, July 1995, |hep-ph/95072^ 
[4] A.Falk and M.Peskin, Phys.Rev.'D49 3320 (1994). 

[5] F.Close, J.G.Kdrner, R.J.N.Phillips and D.J.Summers, J.Phys.GlS 1716 (1992). 

[6] The ALEPH Collaboration, Measurement of the Ab polarization in Z decays, CERN 
PPE/95-156 

[7] G.Bonvicini and L.Randall, Phys.Rev.Lett.73 392 (1994). 

[8] B.Mele and G.Altarelli, Phys.Lett.B299 345 (1993). 

[9] B.Mele, Mod.Phys.Lett.A9 1239 (1994). 
[10] J.K.Kim and Y.G.Kim, Phys.Rev.52 5352 (1995). 
[11] J.Chay, H.Georgi and B.Grinstein, Phys.Lett.B247 399 (1990). 
[12] B.Blok, L.Koyrakh, M.Shifman and A.I.Vainshtein, Phys.Rev.'D49 3356 (1994). 
[13] A.Manohar and M.Wise, Phys.Rev.T>49 1310 (1994). 
[14] J.H.Kiihn, A.Reiter and P.M.Zerwas, Nucl.Phys. B272 560 (1986). 
[15] S.Groote and J.G.Kdrner, MZ-TH/95-17, |hep-ph/950839g 



[16] A.Czarnecki, M.Jezabek, J.G.Korner and J.H.Kiihn, Phys.Rev.Lett.73 384 (1994). 
[17] A.Czarnecki and M.Jezabek, Nucl.Phys.B427 3 (1994). 

[18] A.Falk, M.Luke and M.J. Savage, Hadron spectra for semileptonic decay, UTPT-95-11, 
July 1995. 



22 



[19] P.Ball and V.Braun, Phys.Rev.T>4Q 2472 (1994). 

[20] M.CrisafuUi, V.Gimenez, G.Martinelli and C.T.Sachrajda, First lattice calculation of the 
B meson binding and kinetic energies, CERN-TH- 7521-94, June 1995, |hep-ph/950621C 



[21] Particle Data Group, Phys.Rev.'Db^i 1175 (1994). 
[22] M.B.Voloshin, Phys. Rev. 4934 (1995). 

[23] I.Bigi, M.Shifman, N.G.Uraltsev and A.Vainshtein, Phys.Rev.Tih2 196 (1995). 

[24] I.Bigi, A.G.Grozin, M.Shifman, N.G.Uraltsev and A.Vainshtein, 
Phys.Lett.'BZZQ 160 (1994). 

[25] J.G.Korner and D.Pirjol, Phys.Lett.B334 399 (1994). 



23 



